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High-resolution thermal-expansion and specific-heat measurements were performed on single crys- 
talline BaFe2(Asi_ a! P :E )2 (0 < x < 0.33, x — 1). The observation of clear anomalies allows to 
establish the thermodynamic phase diagram which features a small coexistence region of SDW 
and superconductivity with a steep rise of T c on the underdoped side. Samples that undergo the 
tetragonal-orthorhombic structural transition are detwinned in situ, and the response of the sample 
length to the magneto-structural and superconducting transitions is studied for all three crystal- 
lographic directions. It is shown that a reduction of the magnetic order by superconductivity is 
reflected in all lattice parameters. On the overdoped side, superconductivity affects the lattice pa- 
rameters in much the same way as the SDW on the underdoped side, suggesting an intimate relation 
between the two types of order. Moreover, the uniaxial pressure derivatives of T c are calculated us- 
ing the Ehrenfest relation and are found to be large and anisotropic. A correspondence between 
substitution and uniaxial pressure is established, i.e., uniaxial pressure along the &-axis (c-axis) cor- 
responds to a decrease (increase) of the P content. By studying the electronic contribution to the 
thermal expansion we find evidence for a maximum of the electronic density of states at optimal 
doping. 

PACS numbers: 74.70.Xa, 74.25. Bt, 74.62.Dh, 74.62.Fj 
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I. INTRODUCTION 



In the intensively studied 122 family of iron-based su- 
perconductors, superconductivity can be induced by var- 
ious substitutions in BaFe 2 As 2 (Bal22), e.g., K for Ba, 
Co or Ru for Fe and also P for As, as well as by hydro- 
static pressure^. In all cases, the resulting phase dia- 
grams are surprisingly similar: a superconducting dome 
arises when the magneto-structural transition of the par- 
ent compound is suppressed by either increasing doping 
levels or increasing pressure. The isovalent substitution 
of As by P in BaFe 2 (As 1 _ ;c P :2 .)2 (P-Bal22) with a max- 
imum T c of 30 K is particularly interesting since (nomi- 
nally) no additional charge carriers are introduced. Ob- 
servation of the de Haas-van Alphen (dHvA) effect for 
0.41 < x < 1 (Refs. [7] and [5]) demonstrates the low 
scattering rate of the defects introduced by P substitu- 
tion. Similarly, the absence of quasiparticle scattering 
by the dopant atoms in this system was demonstrated 
in a study on vortex pinning^. It was also found that 
the phase diagram of P-Bal22 can be tracked with hy- 
drostatic pressure from any starting P concentration 10 , 
suggesting that pressure is somehow equivalent to dop- 
ing, as has also been found for the Co-doped systeraHQSl 
Evidence for a quantum critical point close to optimal 
doping was obtained from non- Fermi liquid transport , 
nearly zero Curie- Weiss temperature in NMRpSl and mass 
enhancement in dHvA measurements 7 . 



Here, we present a study of thermodynamic prop- 
erties of the P-Bal22 system focusing on the under- 
doped and slightly overdoped regime (0 < x < 0.33) 
using thermal-expansion, specific-heat and magnetiza- 
tion measurements. Thermal expansion of single crys- 
tals probes the uniaxial pressure dependence of the 
entropy and is therefore well suited to study the ef- 
fect of pressure on superconductivity and spin-density- 
wave (SDW) ordering. Uniaxial pressure effects are ex- 
pected to be quite important in this cla ss of m aterials 
due to their anisotropic crystal structur e 1 12 * 14 * 15 !. High- 
resolution thermal-expansion measurements also provide 
an ideal method for determining the thermodynamic 
phase diagram of these systems. We find that the phase 
diagram of P-Bal22 resembles those of other 122 super- 
conductors, with the exception of a much steeper rise of 
the bulk T c on the underdoped side. We show that the 
small pressure applied by the dilatometer on the sample is 
sufficient to detwin the crystals in situ. By combining the 
results of crystals in the twinned and detwinned states, 
we show that the orthorhombicity of underdoped samples 
is reduced by the onset of superconductivity, which is a 
clear sign of the competition between superconducting 
and SDW/structural order parameters. Moreover, we es- 
tablish a linear relation between P substitution and uni- 
axial pressure, as far as their effects on the phase diagram 
are concerned. Finally, an evaluation of the electronic 
component of the thermal expansion reveals an increas- 
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FIG. 1. Uniaxial thermal-expansion coefficients a.i as a function of temperature T for (a) twinned in-plane and (b) detwinned 
measurements (along the 6-axis) and (c) along the c-axis for various substitution levels as indicated in panel (a). Arrows indicate 
examples of the temperature of the magneto-structural transition T sm and T c ° n . (d) Magnified view of the low-temperature 
region of (a), (e) shows an example of how the splitting of the magneto-structural transition into a structural transition at T s 
and a magnetic transition at Tsdw is inferred from the data. 



ingly pressure sensitive density of states on approaching 
optimal doping from the underdoped side and a maxi- 
mum of the Sommerfeld coefficient at optimal doping. 
Throughout this text, the term "doping" will be used in 
a broad sense, as an equivalent to substitution. 

The paper is organized as follows. In section II and 
III we present the experimental details and results, from 
which we construct the thermodynamic phase diagram 
(section IV). In section V we discuss the interplay of the 
structural/magnetic order with superconductivity. The 
uniaxial pressure effects on T c and the electronic density 
of states are presented in Section VI, and conclusions are 
provided in Section VII. A brief preliminary reporfP^ on 
parts of this work has been presented previously. 



II. EXPERIMENTAL DETAILS 

Single crystals of P-substituted Bal22 were grown 
from stoichiometric mixtures of the starting materi- 
als as described in Ref. [5] Crystals of pure BaFe2P2 
were grown from self- flux using an AI2O3 crucible in 
a closed steel container. Ba and pre-reacted FeP were 
mixed in a ratio of 1:5, heated up to 1300°C and 



slowly cooled down to 1200°C at a rate of 0.3°C/h. 
The typical crystal size is ~ 500 x 500 x 100 /mi 3 . 4- 
circle x-ray structural refinements were conducted on 
three of our samples and yielded a P content of X = 
0.25(1), 0.30(1), 0.33(1). The P content of the other 
samples (x = 0.12(2), 0.18(2), 0.26(2), 0.27(2)) was de- 
termined by energy-dispersive x-ray analysis (EDX) on 
these samples and complemented by 4-circle x-ray diffrac- 
tion on samples from the same batch, hence the larger 
error. 

Thermal expansion was measured in a home-built ca- 
pacitive dilatometer with a typical resolution of 0.1 — 
0.01 A (Ref. [T7)l . In the dilatometer, the sample is 
pressed against one plate of a plate-type capacitor (with a 
force of F « 0.2 N) so that a change of the sample length 
results in a changing capacitor gap. Samples were me- 
chanically detwinned in situ by mounting them such that 
the dilatometer pressure is directed along their tetrago- 
nal [110] direction 1 -^. In this configuration, thermal ex- 
pansion along the (shorter) orthorhombic 5-axis is mea- 
sured. Comparison with twinned samples, for which the 
dilatometer pressure is applied along the tetragonal [100] 
direction, allows to estimate the thermal expansion along 
the orthorhombic a-axis as well. Accurate data could be 
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FIG. 2. The superconducting transition of samples with 
x = 0.27 (underdoped, left column) and x — 0.33 (overdoped, 
right column) seen by different probes. Magnetization (pan- 
els (a) and (b)), specific heat (panels (c)-(d)) and thermal 
expansion (panels (e)-(f)) all show a broad superconducting 
transition with a sharp onset for the underdoped sample and 
more standard, sharp anomalies for the overdoped sample. 



obtained in spite of the (for dilatometry) extremely small 
sample size. 

The specific heat on the tiny samples with masses in 
the order of 100 micrograms was measured with a home- 
made micro-relaxation calorimeter using a 'long relax- 
ation' technique 19 . Each relaxation at different base tem- 
perature provides about 1000 data points over a temper- 
ature interval of up to 50 % above the base temperature. 
Magnetization was measured with a commercial Quan- 
tum Design Vibrating-Sample Superconducting Quan- 
tum Interference Device (VSM-SQUID). 



magneto-structural transition manifests itself as a broad- 
ened peak. Panel (e) shows how this peak can actually 
be decomposed into two separate peaks, revealing a small 
splitting of the structural transition at T s and the mag- 
netic transition at T$dw- We find T s — Tsdw ~ 2.5K for 
x = 0.18 and ~ 7K for x = 0.25 by careful comparison 
of the different in-plane measurements. For x = 0.27 no 
clear splitting of the transition can be observed. Samples 
with x > 0.25 show bulk superconductivity. The onset 
of superconductivity on the underdoped side is signaled 
by a rather sharp kink (at T c on ) followed by a broad 
peak, which has the opposite sign of the anomaly at the 
structural transition. We will use an area-conserving con- 
struction (see panel (e) of Fig. 0) to define a T c mld . On 
the overdoped side (x = 0.33), the anomaly at T c has 
the usual, step-like, shape expected at a second-order 
phase transiton, and T c on and T c mld nearly coincide. 
Note that the sign of the thermal-expansion anomaly al- 
lows us to assign a sample uniquely to lying either on 
the underdoped or on the overdoped side. The anoma- 
lies of the thermal-expansion coefficient along the c-axis 
(panel (c)) have similar shape but opposite sign as com- 
pared to the in-plane measurements in all cases. It is im- 
portant to note that the anomalies of the a^s have the 
same sign and exhibit a very similar shape as repor ted for 
the Ba(Fe 1 _ £C Co x ) 2 As2 (Co-Bal22) systeirPEEU. since 
a,: oc —dS/dpi (uniaxial pressure derivatives of the en- 
tropy), this shows that the entropy of the two systems 
responds similarly to uniaxial pressure. 

The strikingly different shape of the superconduct- 
ing transition of underdoped and overdoped samples is 
also seen in other measurements. Fig. [2] shows the su- 
perconducting transition of the underdoped (x = 0.27) 
and overdoped (a: = 0.33) P-Bal22 samples in the 
dc-magnetization, the electronic specific heat (derived 
by subtracting the lattice heat capacity of an undoped 
Bal22 sample), and the thermal expansion. For under- 
doped BaFe2(Aso.73Po.27)2 there is a sharp onset at T c on 
in all three data sets, however the main anomaly ap- 
pears very broad and rounded. The approximate transi- 
tion width is 10 — 15 K. On the other hand, overdoped 
BaFe2(Aso.67Po.33)2 exhibits sharp, step-like anomalies 
of a, C p and M at T c . A possible origin of the broad- 
ening in the underdoped samples will be presented in 
Section HVl 



III. RESULTS 

Fig. [T] shows the measured uniaxial thermal-expansion 
coefficients a, = (l/Li)dLi/dT where L is the sample 
length and the subscript i stands for the direction. Clear 
anomalies are observed at the magneto-structural tran- 
sition (at T sm ) and at the onset of superconductivity at 
T c on as indicated in the Figure. Panel (b) shows the de- 
twinned measurements of underdoped samples, which be- 
come orthorhombic upon cooling. It is clear that the pro- 
cess of detwinning strongly enhances all anomalies when 
compared to the twinned measurements (panel (a)). The 



IV. PHASE DIAGRAM 

Fig. [3] shows the phase diagram compiled from the 
present thermodynamic data together with previous work 
on resistivity^ of the same system and thermodynamic 
data for Co-Bal22 (Refs. QDandlUD, which have been 
scaled so that the optimal doping concentrations of both 
systems coincide. In underdoped P-Bal22, the bulk T c 
obtained from thermal expansion is lower than T c in- 
ferred from resistivity. Furthermore, the thermodynamic 
measurements reveal a very steep slope of T c . This may 



offer a simple explanation for the broad superconducting 
transitions of the underdoped samples of P-Bal22 (Fig. 
I2J. For a given concentration gradient in the sample, 
the width of the transition will be directly related to this 
slope. For x = 0.25 — 0.27 the width of T c amounts to 
10 — 15 K while it is only ~ IK on the overdoped side. 
This difference can thus be explained by the roughly 13 
times greater slope of T c on the underdoped side. How- 
ever, the sharp kinks at T c on remain unexplained. We 
note that quite broad superconducting transitions on the 
underdoped side have also been reported previously from 
resistivity^ and susceptibility 21 measurements. 

We have already noted that the response of the entropy 
to uniaxial pressure is very similar for P-Bal22 and Co- 
Bal22. Overall their phase diagrams are also very simi- 
lar, however, they also show some important differences. 
First, on the underdoped side, the bulk T c rises more 
steeply in P-Bal22 than in the Co-doped system. This 
may be a signature of a stronger competition between 
the SDW and the superconducting phase. A steep rise of 
the transition line signals that two phases differ little in 
entropy so that only a significant increase of temperature 
triggers the transition. This means that here, supercon- 
ductivity is not much "weaker" than magnetism, and a 
strong effect of the two types of order on each other is 
expected. 

A second difference of the phase diagrams of P-Bal22 
and Co-Bal22 is that the splitting of the magneto- 
structural transition is about two times smaller in the 
P substituted system. Interestingly, the two transi- 
tions remain coincident for isovalent Ru substitution^. 
This suggests that, surprisingly, isovalently substituted 
P-Bal22 is closer to Co-Bal22 than to Ru-Bal22. This 
may be connected to the proposition that a lso Co substi- 
tution does not induce charge doping^l2il The reduced 
magnitude of the splitting may be explained by less scat- 
tering in the P-substituted samples 9 . A reduction of the 
splitting with decreasing disorder has been reported pre- 
viously in the 1111 systems^. Finally, on the overdoped 
side, T c is uniformly higher for P substitution than for 
Co substitution, possibly also due to less disorder. 



INTERPLAY OF ORTHORHOMBICITY, 
SDW AND SUPERCONDUCTIVITY 



Recent studies on Co-Bal22 showed that both mag- 
netic and structural order parameters are reduced in the 
superconducting state^^S. This has been interpreted 
as evidence for microscopic coexistence and competition 
between the orthorhombic SDW phase and superconduc- 
tivity in the underdoped region of the phase diagrarrpSMl 
In some of the 1111 compounds, on the other hand, no 
such coexistence region has been observed^!, and it is 
still debated whether magnetism and superconductivity 
coexist microscopically in K-doped Bal22p™^. it is thus 
of interest to see whether such a coexistence also occurs 
in P-Bal22. A very recent NMR study22 shows evidence 
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FIG. 3. Thermodynamic phase diagrams of P-Bal22 (lower 
axis) compiled from the present thermal-expansion data and 
of Co-Bal22 (upper axis) scaled so that the maximum T c 's co- 
incide. Data for Co-Bal22 was obtained from specific-hea1p3 
and thermal-expansion 12 measurements. Also added are pre- 
viously published resistivity data 5 on the P-Bal22 system. 
Clearly, the bulk T c rises more steeply for P-Bal22 than for 
Co-Bal22 on the underdoped side. 



for a reduced magnetic order parameter in the super- 
conducting state also for P-Bal22 and therefore claims 
microscopic coexistence of magnetism and superconduc- 
tivity. In the following we will use our thermal-expansion 
data to study this interplay between structural/ magnetic 
order and superconductivity via the response of the sam- 
ples' lengths and the orthorhombic order parameter to 
superconductivity. 

As an example, Fig. [4] shows the relative length 
changes of the x — 0.25 sample. We note that the a- 
axis length was not directly measured. However, it can 
be derived from the difference between in-plane twinned 
and detwinned data. The difference of the relative length 
changes of a- and 6-axes, i.e., the orthorhombic order pa- 
rameter S(T) inferred from our measurements, is shown 
in the inset of Fig. |4j We see that a- and &-axis 
lengths start to differ clearly at the structural transi- 
tion, and this difference is directly proportional to the 
structural order parameter in this system. There exists 
a high-temperature tail to the transition, which arises 
from the small applied in-plane pressure of the order of 
5 — lOMPaPSBH Clearly, a- and 6-axes approach each 
other again below T c and the orthorhombic order param- 
eter decreases, suggesting a similar coupling of structural 
and superconducting order parameters as observed in Co- 
Bal22 (Ref. 28). Interestingly, the c-axis length increases 
below T s and then decreases below T c , and thus exhibits 
a very similar behavior as the in-plane axes. Although it 
is not clear how, or if at all, the c-axis response can be 
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FIG. 4. Relative change of sample length of the 

BaFe2(Aso.75Po.25)2 sample, measured in the twinned (in- 
plane average, black line) and detwinned (b-axis, red line) 
configuration and along the c-axis (blue line). The a-axis 
length (broken purple line) was estimated from the twinned 
and detwinned in-plane measurements and the volume change 
(green broken line) from the twinned in-plane and the c-axis 
measurements. The inset shows the difference of a-axis and 
6-axis. For this sample T sm « 55 K and T c mld = 16 K. 
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FIG. 5. Relative change of the 6-axis length (in-plane length 
for the overdoped x = 0.33 sample) for various doping levels. 
Its shortening signals the structural transition and its increase 
below T c for the superconducting samples demonstrates the 
competition between the orthorhombic phase and supercon- 
ductivity. The upper inset shows the data for x = 0.33 on 
an expanded vertical scale. The lower inset shows the doping 
evolution of dbtot as defined in the main panel. 



linked to the structural/magnetic order parameters, our 
results directly show that the SDW state favors a longer 
c-axis. Equivalently, compressing the c-axis by uniaxial 
pressure will destabilize magnetism. This is in agree- 
ment with DFT calculations under uniaxial pressure^. 
The decrease of the c-axis length upon entering the su- 
perconducting state is thus consistent with a suppression 
of magnetism by superconductivity. Finally, we note that 
the effect of these transitions on the volume is very small 
(see Fig. B| due to an almost complete cancellation of 
the anomalies along the different directions. 

In Figures [5] and [6] we show the length changes as a 
function of temperature for different doping levels for 
the 6- and c-axes, respectively. These axes can be mea- 
sured accurately. We do not show the twinned and in- 
ferred a-axis data or 6(T), since these are not as reliable 
due to possible partial detwinning during the twinned 
in-plane measurements. We note that the tempera- 
ture dependence of the orthorhombic order parameter 
S = (a — b)/(a + 6) is however reflected in the (neg- 
ative) 6-axis length since a- and b- axis are found to 
evolve nearly symmetrically. In Fig. [5] we define a pa- 
rameter Abtot, which quantifies the total change of the 
6-axis length due to the magneto-structural transition. 
As shown in the inset of this figure, Ab to t decreases lin- 
early with P content confirming reliable detwinning. An 
exception are the samples with x — 0.26 and 0.27. Even 
though the transition temperatures are similar, Abtot is 
much smaller for these samples and does not follow T s 
any more. Possibly, only a part of the sample undergoes 
the magneto-structural transition. Similarly, NMR found 



evidence that a fraction of a collection of samples with 
x = 0.25 does not undergo the SDW transition 21 . Hence 
we exclude the samples with x — 0.26 and 0.27 in the 
analysis of the electronic part of the thermal expansion 
(Section |VI[ ). As already seen for the x = 0.25 sample, the 
6-axis length increases and the orthorhombicity is signif- 
icantly reduced when samples become superconducting. 
On the other hand, the in-plane length of the overdoped 
x = 0.33 sample decreases below T c . However, the effect 
is tiny in comparison to the underdoped samples (see in- 
set in the upper left corner of Fig. IHJ. 

Interestingly, the c-axis length displays a similar be- 
havior as the in-plane axes. Fig. [6] (a) shows length 
changes of the c-axis in analogy to Fig. [5] The iron- 
based materials are more compressible along the c-axis 
than in-plane, which is reflected in a larger lattice con- 
tribution to the thermal expansion along the c-axis^HIl 
To obtain the "electronic c-axis length changes" ALf a 
suitable lattice background has to be subtracted. For 
this lattice background we used the x — 1 data, because 
pure BaFe2P2 does not undergo any phase transitions 
and, additionally, the electronic contribution to its ther- 
mal expansion is negligibly small compared to the ef- 
fects discussed here. We assume that this background 
is independent of P doping and, furthermore, that the 
high-temperature value of the total cv c is doping inde- 
pendent. Variations of this high-temperature value are 
attributed to uncertainties arising from the tiny sample 
lengths along c (~ 70 — 150 fim). These deviations can 
be modeled with a correction factor between 0.9 and 1.3 
to the lattice contribution for the different samples. This 




80 

T[K] 



40 80 120 160 

T[K] 



FIG. 6. (a) Relative change of the c-axis length for various 
doping levels. Panel (b) shows the data after subtraction of 
a background (the data for x — 1 times an individual factor 
close to 1, see text). 



correction does not affect the shape of the anomalies. 
The thus obtained ALf are shown in panel (b) of Figure 
M as a function of temperature. Intriguingly, these ALf 
resemble closely the T dependence of the orthorhombic 
order parameter for all doping levels. This opens the 
interesting possibility to study the orthorhombic SDW 
state via its c-axis length. As already detailed in the 
analysis of the x = 0.25 sample, it is clearly seen that 
the magnetic state favors a longer c-axis. Its destabi- 
lization by superconductivity is reflected in a reduced 
c-axis length. Interestingly, along the c-axis, there is a 
marked response also for the overdoped x = 0.33 sample. 
The superconducting state in the absence of static mag- 
netic order apparently favors a longer c-axis in a similar 
manner as does the SDW state of underdoped P-Bal22. 
The same effect, though much smaller in magnitude, is 
also seen for the in-plane dimension: The a-axis length 
of overdoped P-Bal22 is (slightly) reduced in response to 
superconductivity, as is the average in-plane length of the 
underdoped material upon entering the SDW state. The 
similar response of the lattice parameters may hint at a 
kinship of SDW and superconductivity as is suggested by 
nearly identical Griineisen-parameters of the two phases 
in Co-Bal22SH. 

Finally, we address the question of microscopic coexis- 
tence of magnetism and superconductivity. Dilatometry 
is a macroscopic probe and the T dependence of a sample 
length cannot a priori be equated with the T dependence 
of the lattice parameters if one has a phase-separated 
sample. Still, the close resemblance between our results 
and the evolution of lattice parameters of Co-Bal22 ob- 
served by x-ray diffraction 2 ^ supports that such an iden- 
tification may indeed be done, suggesting that our data 
really show the effect of superconductivity on the lat- 
tice constants. This effect is most naturally explained by 
homogeneous coexistence and competition of magnetism 
and superconductivity also in P-Bal22. In confirmation, 
NMR measurements 21 show an influence of superconduc- 



tivity on the local magnetic order parameter and hence 
also conclude that SDW and superconductivity coexist 
microscopically but compete with each other. 

Summarizing this section, the temperature and doping 
dependence of the lattice parameters of underdoped P- 
Bal22 reflect how superconductivity competes with the 
orthorhombic magnetic state and finally suppresses it. 
In contrast, superconductivity on the overdoped side is 
stabilized by lattice changes similar to those that occur 
upon entering the SDW phase, which hints at a close 
similarity of the two ordering phenomena. 



VI. UNIAXIAL PRESSURE EFFECTS 

For second-order phase transitions, uniaxial pressure 
derivatives of the transition temperature, dT c /api, can be 
deduced from the jumps of the thermal-expansion co- 
efficient Acq and the specific heat AC p . For the over- 
doped, tetragonal, sample (x = 0.33, T c — 29 K), clear 
second-order jumps in the uniaxial thermal-expansion co- 
efficients Aa a = 0.9(2) x 10- 6 K-\ Aa c = -4.0(5) x 
lO^K -1 and the specific heat AC P /T C = 38(4) mJ mol -1 
K~ 2 are observed (see Fig. 2b. The magnitude of the 
specific-heat jump is in good agreement with data of a 
recent report by Chaparro et al! 3 -^. dT c /dpi can be calcu- 
lated via the Ehrenfest relation 



dpi 



V„, 



Mi 



AC P /T C 



(1) 



where V m = 59.6cm 3 /mol (Ref. [S} is the molar volume. 
The data yield aT c /ap a = 1.4(5) K/GPa and aT c /dp c = 
—6.3(1.2) K/GPa. The hydrostatic pressure derivative 
dT c /dp = 2aT c /dp a + aT c /dp c = -3.5(2.2) K/GPa is in 
reasonable agreement with high-pressure experiments 38 
which report an initial slope of oT c /dp — —1.8 K/GPa for 
a sample with T C =30.5K. The same kind of anisotropy 
and a similar magnitude of the uniaxial pressure deriva- 
tives have been observed for a slightly overdoped Co- 
Bal22 crystaP^. The pressure derivatives on the un- 
derdoped side are more difficult to extract due to the 
large widths of the transitions. However, it is clear that 
the dTc/dpi's are of opposite sign and by far larger in 
magnitude than those of the overdoped sample. Using 
the construction for the jumps shown in Fig. 2(e), we 
find for the x = 0.27 sample dT c /dp ab = -18(8) K/GPa, 
aT c /dp b = -44(17) K/GPa and aT c /dp c = 20(8) K/GPa. 
Here, p ab refers to an average of pressure along the a- 
and 5-axis, as is relevant to the twinned in-plane mea- 
surements. 

Since AC P is always positive, the sign of the anomalies 
in the a* provide information on the anisotropy of the 
dT c /dpi even in the absence of specific-heat data. Strik- 
ingly, the uniaxial pressure derivatives of T s and T$dw 
have the opposite sign to the derivatives of T c for all 
directions^. This is again consistent with a competi- 
tion between the orthorhombic SDW phase and super- 
conductivity, since pressure will favor either the SDW or 
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FIG. 7. The electronic contribution to the uniaxial thermal-expansion coefficients divided by T, a e i/T, obtained from sub- 
tracting the data of the x — 1 sample as a phonon background (a) in the twinned in-plane measurements, (b) along the 
orthorhombic 6-axis (detwinned in-plane measurements) and (c) along the c-axis. Arrows indicate how values of the normal 



state af/T(T — 0), which are reported in Fig. ^l are extracted. 



superconducting phase at the expense of the other. In- 
terestingly, all signs of the derivatives can be accounted 
for by identifying uniaxial pressure with a shift in the P 
content: stress along the c-axis (and a-axis on the un- 
derdoped side) corresponds to an increased P content 
while stress along the 6-axis (averaged in-plane axis) cor- 
responds to a lower P content. This means that the phase 
diagram may be tracked (forwards or backwards) by the 
application of uniaxial pressure just as it can be tracked 
by hydrostatic pressure^. By comparing the values of 
dT c /dpi obtained above with the respective slopes in the 
phase diagram dT c /dc — 6.7K/at%(P) (underdoped) and 
dT c /dx = — 0.5K/at%(P) (overdoped) one can estimate 
that increasing the P content by 1 at% corresponds to 
a uniaxial pressure of — 0.35GPa (— 0.15GPa) in the in- 
plane average (along the 6-axis). The lower accuracy of 
the c-axis data unfortunately does not allow such a quan- 
titative analysis. 

Thermal expansion data can also be used to obtain 
information about the (uniaxial) pressure dependence 
of the electronic density of states. For a Fermi liquid, 
the electronic thermal expansion divided by T, af/T, 
is expected to be a constant, a e i t i/T = — (l/V m )dy/dpi, 
with 7 the Sommerfeld coefficient. (This follows directly 
from cti — —(l/V m )dS/dpi and the electronic specific heat 
C — 7T.) Fig. [7] shows the electronic/magnetic thermal- 
expansion coefficients af/T ex —dy/dpi, which were ob- 
tained from the original data by subtracting the data 
for x = 1 as a phonon background using no correction 
factor (see section Iv}. With increasing doping, a sizeable 
contribution from electronic/magnetic degrees of freedom 
evolves below the SDW transition for all three crystal 
directions. Especially for the x = 0.18 data, this contri- 
bution is nearly independent of temperature up to about 



60 K, consistent with Fermi-liquid-type behavior. This is 
reminiscent of what is observed in weak itinerant ferro- 
magnetic systems, such as MnSi (Ref. HCT|) . and points 
to an intricate mixing of highly pressure-dependent elec- 
tronic and magnetic degrees of freedom. For x = 0.25, 
Tsdw & n d T c are too close to each other so that this 
constant term in a e i/T is not seen. 

Panel (a) of Fig. |] shows af/T(T = 0), obtained 
by extrapolating the normal state af/T data to zero, 
and corresponding dy/dpi values versus P substitution. 
There is a sharp increase in the magnitude of dy/dpi 
upon approaching optimal doping and a sign change 
when passing to the overdoped side. The discontinuity 
of the derivatives dy/dpi at optimal doping is also sug- 
gested from an analogy with the Co-doped systenfSESl 
We can obtain 7 as a function of x from integration of 
these data if, as before, we equate doping with pressure. 
We recall that we can link average in-plane (c-axis) pres- 
sure to a decreased (increased) P content. As is evi- 
dent from the sign change then, 7(3;) has a maximum 
at optimal doping. More quantitatively, the proportion- 
ality factor from above (dp a b/dx = — 0.35GPa/at%(P) 
and dpb/dc = — 0.15GPa/at%(P)) yields values for dy/dc 
as a function of doping; for example dy/dx = 2.4 mJ 
mor 1 K~ 2 /at%(P) at X = 0.25. Integration of the 
smooth lines in Fig. [8ja) rescaled by this factor results 
in j(x) shown in Fig. p[b). 7 increases with doping 
in the underdoped region, has a sharp maximum at opti- 
mal doping and then decreases with further doping. Also 
shown are measured valuetP^of 7, which agree reasonably 
well considering the assumptions and uncertainties of our 
approach and the uncertainties in the values of 7. The 
overall behavior of 7(2:) is very similar to what has been 
observed in Co-doped Bal22, where the initial increase 
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FIG. 8. (a) Values of the electronic thermal expansivity 
ttel,i/T extrapolated to zero temperature. The right-hand 
scale shows the corresponding values of cty/dpi. Lines are a 
guide to the eye. (b) The Sommerfeld coefficient 7 obtained 
by integrating the smooth lines from (a) (black line), see text 
for details. Full squares indicate directly measured values of 
7 and open squares indicate values estimated from the mea- 
sured specific-heat jump^H, see also note I39I 



VII. CONCLUSIONS 



We have constructed the thermodynamic phase dia- 
gram of P-Bal22 in a significant range of P substitu- 
tion < x < 0.33. Compared to the Co-substituted 
system, the superconducting T c rises more steeply and 
there is only a small overlap region where superconduc- 
tivity occurs within the orthorhombic SDW phase. P- 
Bal22 seems to be rather close to a first-order transition 
between magnetism and superconductivity, as is realized 
in some of the 1111 systems. The question arises why 
P-Bal22 and Co-Bal22 differ in this aspect of interplay 
between SDW and superconductivity. The competition 
of the two types of order is also reflected in the lattice 
constants. The destabilization of the SDW by super- 
conductivity is reflected both in the in-plane lengths, as 
a reduction of orthorhombicity, and in the c-axis length. 
However, superconductivity on its own favors a change of 
lattice parameters similar to the SDW state, suggesting 
an intimate relationship of the two ordering phenomena. 
Uniaxial pressure derivatives of the transition tempera- 
tures are large and anisotropic. Considering the effects on 
the phase diagram, we put forward a linear relation be- 
tween uniaxial pressure and P content. The Sommerfeld 
coefficient is increasingly pressure (and therefore doping) 
dependent on approaching optimal doping. By evaluat- 
ing its derivatives, we deduce that the density of states 
has a maximum at optimal doping. 
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